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Abstract. - We show that the inclusion of effects beyond Migdal's limit in the electron-phonon 
interaction naturally leads to an isotope effect for the effective mass m* of the charge carriers 
even much before reaching the polaron limit. This is the situation already considered in our 
approach to nonadiabatic superconductivity [j^]. Such a result provides a scenario different 
from the polaronic one for the interpretation of the recently observed isotope effect on m* in 
YBa2Cu306+a: and La 2 -^Sr :E Cu04 0, f|. 



The recent observation of an oxygen-mass dependent penetration depth A(0) in YBa2Cu3 06+: 
[0 and in L^-zSr^CuO^ j|] has raised the question whether this effect could be attributed 
to the breakdown of the Migdal approximation [||. In fact, according to the classical Migdal- 
Eliashberg theory of superconductivity, the zero temperature penetration depth in the London 
limit and in the absence of magnetic impurities Q is given by: 



A(0)oc A / — , (1) 

V ^ 

where to* is the electronic effective mass and n s is the supercarrier density. According to eq.(Q), 
an isotope effect on A(0) can therefore be induced by an isotope effect on to* and/or on n s . 
However, for T <C T c , also assuming an electron-phonon (el-ph) type of pairing, the dependence 
of n s on the ion- mass M is negligible |^| and a possible isotope-induced dependence of n s on the 
hole concentration has been ruled out 0|. Zhao et al. therefore concluded that the observed 
isotope effect on A(0) can be entirely attributed to an oxygen- mass dependent m*. In this 
hypothesis and by using eq.(Q), the oxygen isotope effect on m*, a m - = —dlii(m*)/d\n(M), 
has been estimated to be a TO » = —0.61 ±0.09 for YBa 2 Cu3 06.94 0] and a m * ~ —0.8 (—0.5) for 
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La 2 _a;Sr x Cu04 for x = 0.105 (0.15) ||]. These large negative values of a m * cannot be explained 
by the classical Migdal-Eliashberg (ME) theory of the el-ph interaction, which states that 
ot m * — 0. In fact this theory predicts that the electron-mass renormalization factor m* /m = Zq 
is given by Zq = 1+A, where A is the el-ph coupling constant, i.e., a quantity independent of the 
ion-mass M (A must not be confused with A(0), the penetration depth). This is a consequence 
of the adiabatic condition Xujjj/igvp) <C 1, where lojj is the Debye phonon frequency, Vp is the 
Fermi velocity and q is the momentum transfer in the el-ph scattering. The adiabatic condition 
permits to neglect the vertex corrections to the el-ph interaction (Migdal approximation) and 
represents the basis of the formulation of the ME theory. The unusual isotope effect a m * 
observed in YBa2Cu306+a; and in La2- :E Sr :E Cu04 could therefore be interpreted as a clear 
signature of the breakdown of the Migdal approximation. 

However, from the quantitative point of view, some care has to be used in estimating 
a. m * from the isotope effect of A(0) via eq.(|l|), since the latter applies only when the Migdal 
approximation is a valid assumption. A generalization of eq.(^) beyond Migdal's theorem 
requires in fact the knowledge of the Eliashberg equations for temperatures below T c with 
the inclusion of the nonadiabatic contributions, i.e., a task which even in the perturbative 
approach presents serious technical difficulties. Nevertheless, the measurements reported in 
Refs.|j], ||] lead to two important consequences. First, the observation of an isotope effect on 
A(0) implies a significant role of the el-ph interaction in the high-T c compounds. Second, as 
pointed out before, such an the el-ph interaction falls outside of the validity of the Migdal 
approximation. Therefore, the findings of Refs.jl], ^) open new perspectives in investigating 
unusual isotope effects on quantities which, according to the classical (in the sense of ME) 
theory of the el-ph interaction, should not show isotope effects at all. 

One of the most unambiguous evidences for the breakdown of the ME theory could be 
provided by the observation in the normal state of an isotope effect on the effective electronic 
mass to*. In principle, isotope sensitive specific heat measurements could be able to observe 
such an effect provided that the electronic contribution to the specific heat can be clearly 
singled out ||. With this motivation, we consider in this paper the consequences of the 
breakdown of Migdal's theorem on m* and estimate the coefficient a m * of the electron-mass 
isotope effect. We show that nonzero values of a m * in the normal state can be obtained by 
considering two different regimes of the el-ph coupled system: the polaronic state and the 
"nonadiabatic regime" as described by the inclusion of the first el-ph vertex correction in the 
electronic self-energy. We have already studied the latter situation in connection with the 
superconductive transition ]?], || and here we show that this theory, which does not imply the 
crossover towards the polaronic state, naturally leads to an isotope effect for m*. 

Let us start our discussion with the polaron model. A similar analysis has already been 
performed in Ref. j^j in connection with the experimentally observed isotope effect on the 
penetration depth A(0) {{],§]. For simplicity, we consider the Holstein model [^) where phonons 
with frequency loq are locally coupled to electrons through a structureless el-ph coupling g 
(small polaron). A polaronic state is characterized by a strong electron-lattice correlation 
[ELO] and/or a large effective mass @. Within the range of applicability of the Holstein 
approximation, the effective polaron mass is given as follows: 



m* ~ mexp ( —n ) , (2) 

where m is the bare electron mass. Equation (^|) applies to the el-ph system in the antiadiabatic 
regime loq/Ef > 1, provided that the quantity g 2 /wo is not too large [[H]. Since g 2 oc (Mluq)^ 1 
and Wq oc 1/M, from eq.(S) we obtain that to* oc exp \J^M, where 7 is a constant independent 
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Fig. 1. - Electronic self-energy beyond Migdal's approximation. Solid lines represent electronic 
Green's functions while the wavy lines are phonon propagators. 



of M. The isotope effect on m* is therefore given by: 



dln(m*) 1 / g 



d\n(M) 



2 ' 

(3) 



The above expression gives a negative value of a m * in accordance therefore with the results 
reported in Rcfs.|l], |2). We must stress however that this result is based on eq.(||), which 
holds true only as long as u>o/Ef > 1, i. e., a rather inadeguate limit for Yl^CuaOe+z 
and La2- a: Sr a ;Cu04 |L2|. It is interesting to notice that, at optimal doping, LSCO shows 
a large negative a m * [0 and a negligible isotope coefficient of the critical temperature ar c 
[ ft3| . This behavior is therefore in contrast with the prediction of the bi-polaronic theory of 
superconductivity which claims that a m * oc ax c ■ 

Besides the ME and polaronic scenarios, the el-ph coupled system may display also a regime 
which is beyond Migdal's limit (and therefore beyond ME framework) and well separated 
from the crossover between the quasi free-electron and the polaronic state. Such a regime, 
which we call nonadiabatic, is characterized by quasi free-electron states (A <^ 1) coupled 
in a nonadiabatic way to the lattice so that the el-ph vertex corrections are relevant. In 
the nonadiabatic regime, the electronic self-energy Y,(ioj n ) is given by the graph depicted in 
Fig. 1 where the first vertex correction has been included. The mass renormalization factor 
Zq = m*/m is obtained by the iuj n — > limit of Z(iuj n ) = 1 — Y>(iui n ) / (iuj n ) . The function 
Z(iu) n ) has been obtained in a previous work and its expression is reported below: 



Z(tuj n ) = 1 H > — r - arctan — — . (4) 



Here, uj n and uo m are fermionic Matsubara frequencies and A^(iw„, ico m ; Q c ) is the frequency 
dependent el-ph coupling resulting from the inclusion of the first vertex correction function 
Pv{itu n ,iuj m ;Q c ) into the electronic self-energy: 



X Z (i0J n , i^m] Qc) = A [1 + XP V {iLOn, i0J,n, Qc)} ■ (5) 

The explicit expression of the vertex function Py has already been presented in Ref.|| for 
the three dimensional case and in Ref. Jl4]| for the two dimensional one. The dimensionless 
parameter Q c — q c /(2kp), where q c is a cut-off over the momentum transfer q and Uf 
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Fig. 2. - Renormalization function Z(iu) for luq/Ef — 0.2 and A = 1. Solid lines: case with 
vertex correction for Q c = 0.1, 0.2, . . . , 0.5. Dashed lines: case without vertex correction. 

is the Fermi wave number, follows from the model we use for the el-ph coupling function 
9 2 (Q) — (q 2 /Qc)0(Qc — Q), where Q = q/(2kp) and 9 is the Heaviside function. The aim 
of this model is to simulate the effect of strong electronic correlations on the el-ph matrix 
element. In fact, according to different theoretical approaches, the tendency of the electronic 
correlation |L5| and the weak screening Jl6[ is to suppress the scattering processes with large 
momentum transfer and, at the same time, to enhance the small Q — q/(2kp) scatterings. 

The calculated Z(iu>) is shown in Fig. 2 as a function of ui/uiq for uoq/Ef — 0.2, A = 1.0. 
When the vertex correction is not included (dashed line), Z(iu) takes the maximum value 
Zmax at u> = 0. The deviation from the adiabatic limit [lim^ /e f —>o Zmax = 1 + A] is due 
only to the finiteness of the electronic band ||. When the vertex correction is included 
(solid lines) the maximum of Z(iuj) is shifted towards higher values of the frequency u). For 
fixed values of loq/Ef and A, the position and the amplitude of Z max depend on the cut-off 
parameter Q c . Moreover, the value of Z{ilo) at u> — is considerably lowered by the presence 
of the vertex correction. This feature can be understood by considering that at u> — the 
electron-mass renormalization factor Zq is mostly modified by the static limit of the vertex 
function, Py{iuj n ,iuj m — * iuj n ;Q c ) , which is found to be negative [^J. This can be recovered 
also by the following analytic expression valid for uiq/Ef — > and cuq/ '(Q 2 c Ef) small: 

where the contribution of the negative static limit is given by the third term of the right hand 
side. Note instead that for the superconducting transition the situation is rather different, 
since the range of the relevant frequencies is of order of ojq. In such a region of frequencies the 
vertex function shows a complex behavior that can lead to an amplification or a suppression of 
T c depending on the value of the parameter Q c R ^ . This different role of the nonadiabatic 
contribution in Zq and T c reflects the fact that the inclusion of the vertex corrections cannot 
be merely interpreted as a simple renormalization of the el-ph coupling. 

We show in Fig. 3 the calculated coefficient of (Xtyi* ELS cL function of the adiabatic parameter 
ujq/Ef for A = 1 and different values of Q c . a m * takes negative values regardless the absence 
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Fig. 3. - Coefficient a m * of the isotope effect of the effective electronic mass calculated for 
A = 1. Solid lines: case with vertex correction for Q c = 0.2, 0.3, . . . , 0.6. Dashed lines: case 
without vertex correction. 



(dashed line) or the presence (solid lines) of the vertex correction. The presence of a nonzero 
a m » for the case without vertex correction has to be ascribed only to finite-band effects, 
which give rather small absolute values of the isotope coefficient. The inclusion of the vertex 
correction amplifies the nonadiabatic effects and leads to more negative values of a m * with 
respect to the case without vertex correction. This can be also inferred from eq.(^|) which 
gives: 



2Z 4 \Q 2 C E F 



w 



(7) 



where the finite band effects have been neglected {ujq/Ef <C 1) and the vertex correction is 
evaluated only up to the linear term ujq/ '(Q^Ep). Note that in Fig. 3 a m - does not show 
in the whole region of loq/Ef a crucial dependence on Q c as instead is observed for the 
superconducting transition temperature T c [jjj |j. As a consequence, the isotope effects on T c 
and m* are not proportional as in the bi-polaron theory of superconductivity and in principle 
it is possible to have a negligible ax c and at the same time an appreciably nonzero a m - . 

The result reported in Fig. 3 have been obtained by considering a structureless electronic 
density of states (DOS) and are therefore relevant for three dimensional systems like the 
fullcrene and the BaBiC>3 compounds. However, photoemission measurements have given 
evidence for flat electronic bands close to the Fermi level for a large class of cuprates [ p"7[ , 
opening the possibility that van Hove singularities (vHs) in the DOS could lead to important 
effects on both the normal and the superconducting phases []l8| . As discussed in Ref. [|l4| , the 
presence of a vHs near the Fermi level represents an intrinsic nonadiabatic situation, where 
Migdal's theorem is not valid even if loq/Ef « 1. It is natural therefore to investigate the 
effect on a m * of the vertex correction calculated by using a vHs in the DOS. Here we anticipate 
some preliminary results valid for uq/Ef <C 1 obtained by using the following DOS: 



N(e) = -No In 



(8) 
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where —E/2 < e < +E/2 and iVo = N/E corresponds to a constant DOS with N electronic 
states. In this model, the vHs is located at the Fermi energy Ep = E/2. By making use of 
cq.(^), we have found that for Ao = 1 and ujq/Ef = 0.05, the isotope coefficient a m * is —0.29 
(—0.6) for Q c = 0.4 (0.2). Therefore, the combined effect of the vHs and the vertex correction 
leads to a magnification of the values of a m * , in better agreement with the experimental ones. 

In summary, we have discussed the effect of different el-ph regimes on the effective electronic 
mass to* . The main result is that the possibility of having nonzero values of a m * cannot be 
associated exclusively to the presence of polaronic charge carriers. We have shown, in fact, 
that a negative ro*-isotope coefficient can be obtained by taking into account the first el-ph 
vertex correction beyond Migdal's limit. The latter result appear to be of particular interest 
in view of the fact that in the nonadiabatic theory of superconductivity the breakdown of 
Migdal's theorem can lead to a strong enhancement of T c and various other effects (7|, [|, Q . 
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REFERENCES 



Zhao G. M. and Morris D. E., Phys. Rev. B, 51 (1995) 16487. 
Zhao G. M., Hunt M. B., Keller H. and Muller K. A., Nature, 385 (1997) 236. 

Migdal A. B., Sov. Phys. JETP, 7 (1958 ) 996. 

Bill. A, Kresin V. Z. and Wolf S. A., |cond-mat/9801186| Preprint (1998). 
Abrikosov A. A., Gorkov L. P. and Dzyaloshinski I. E., Methods of Quantum Field Theory 
in Statistical Physics, (Prentice-Hall, Englewood Cliffs) 1963. 
Loram J. W. et al: ., Physica C, 235-240 (1994) 134. 

Grimaldi C, Pietronero L. and Strassler S., Phys. Rev. Lett., 75 (1995) 1158. 
Pietronero L., Strassler S. and Grimaldi C, Phys. Rev. B, 52 (1995) 1995; Grimaldi C, 
Pietronero L. and Strassler S., Phys. Rev. B, 52 (1995) 10530. 
Holstein T., Ann. Phys., 8 (1959) 343. 

Capone M., Stephan W. and Grilli M., Phys. Rev. B, 56 (1997) 4484. 
Ciuchi S., de Pasquale F., Fratini S. and Feinberg D., Phys. Rev. B, 56 (1997) 4494. 
Plakida N. M., High- Temperature Superconductivity,, (Springer, Berlin) 1995. 
Franck J. P. et al: , Phys. Rev. Lett, 71 (1991) 283. 

Cappelluti E. and Pietronero L., Phys. Rev. B, 53 (1996) 932; ibid., Europhys. Lett., 36 
(1996) 619. 

Zeyher R. and Kulic M., Phys. Rev. B, 53 (1996) 2850; Grilli M. and Castellani C, Phys. 
Rev. B, 50 (1994) 16880. 

Weger M., Peter M. and Pitaevskii L. P., Z. Phys. B, 101 (1996) 573. 

For an extensive review see: Shen Z. X. and Dessau D. S., Phys. Rep., 253 (1995) 1, and 
references therein. 

[18] Labbe J. and Bok J., Europhys. Lett., 3 (1987) 1225; Newns D. M. et al: ., Comments 
Condens. Matter Phys., 15 (1992) 273; Abrikosov A. A., Physica C, 244 (1995) 243. 



